Compact surface Fano states embedded in the continuum of waveguide arrays 
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Wave localization at surfaces is a fundamental 
physical phenomenon with many important ex- 
amples such as surface plasmon polaritons [1, 2j, 
Dyakonov surface modes [3j , and Tamm states [4 J . 
Although being of different origin, these modes 
share the same property: they decay exponen- 
tially away from the surface, with the eigenfre- 
quencies of the surface states being outside the 
band of the spatially extended propagating states. 
Here we describe theoretically and observe ex- 
perimentally a novel type of surface states in a 
semi-infinite waveguide array with simultaneous 
Fano and Fabry-Perot resonances. We demon- 
strate that the Fano surface mode is compact, 
with all energy concentrated in a few waveguides 
at the edge and no field penetration beyond the 
side-coupled waveguide. We show that by broad- 
ening the spectral band in the rest of the waveg- 
uide array it is possible to suppress exponentially 
localized modes, whereas the Fano state having 
the eigenvalue embedded in the continuum is pre- 
served. 

Remarkably, in their seminal paper von Neumann and 
Wigner [5] explicitly constructed a potential that sup- 
ports so-called "bound states in the continuum" (BIC), 
i.e., a particular localized mode with the eigenfrequency 
inside the band of extended states. Since then, the pecu- 
liar concept of BIC attracted a lot of attention in various 
branches of physics [6j 0, reaching a climax with the 
direct observation of an electronic bound state above a 
potential well due to Bragg reflection in semiconductor 
heterostructures [8]. In the optical domain it was shown 
that BIC can be generated in photonic crystals and opti- 
cal waveguide arrays [9-16 . The first true observation of 
a BIC in any system was carried out in an optical system: 
a planar optical waveguide array with two side-coupled 
defects, in which an in-band defect state is formed by 
virtue of symmetry [17]. Recently, a special design of a 
one-dimensional (ID) waveguide array with modulated 
inter-site coupling was suggested in order to obtain "sur- 
face BIC" HBJ. 

The phenomenon of Fano resonance [19] provides a 
simple and universal means for obtaining bound states in 
the continuum, which can be realized in optics through 
the addition of side-coupled waveguides or cavities [20] . 
Indeed, at Fano resonance there appears an effectively 
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FIG. 1. (Color online) Sketch of a semi-infinite waveguide 
array, where a side- coupled waveguide is introduced to provide 
the Fano resonance. 



infinite potential barrier such that an incident wave ex- 
hibits complete reflection, and such mechanism is re- 
sponsible for light trapping studied in photonic crystal 
waveguides between two side coupled defects [2TH23] , at 
a cavity placed near the waveguide termination [24 , and 
predicted for waveguide arrays with side-coupled wave- 
guides [9 . In this regime the mode can be trapped even 
when its eigenvalue is embedded in a continuum spec- 
trum of the optical waveguide. 

In this Letter, we present the theoretical analysis and 
the experimental observation of the formation of surface 
Fano states. To this end, we employ an optical system: 
an array of evanescently coupled waveguides with a side 
coupled waveguide close to the lattice edge. We show 
that the Fano modes are "compact", i.e., there is no 
wave penetration beyond the location of the side-coupled 
waveguide, in agreement with the principles of Fano 
resonance based on complete destructive interference in 
transmission. Based on this property, we demonstrate 
that usual exponentially localized defect modes can be 
suppressed by increasing the waveguide coupling past the 
side-coupled waveguide: These modes disappear when 
their eigenvalue enters the broadened spectral band, in 
contrast to the Fano mode that remains trapped. Such a 
selective manipulation of modes sensitive to their local- 
ization mechanism may find applications in various phys- 
ical contexts, since Fano resonance is a general wave phe- 
nomenon. We study the formation of Fano surface local- 
ized modes in a planar optical waveguide array with one 
side-coupled waveguide as shown in Fig. [I] In analogy to 
the atomic physics context where Fano results were orig- 
inally derived, this side-coupled waveguide is called the 
"autoionizing" waveguide (AW). The waveguide struc- 
ture consists of a ID chain of identical waveguides. The 
propagation direction and the transverse direction are z 



and x, respectively. The chain is coupled to the AW 
with strength V. The part containing the waveguides 1 
to N is a resonator region. Here, the coupling coefficient 
between neighboring waveguides is C. We consider the 
general situation when the coupling strength Cb in the 
"background" region, i.e., the semi- infinite array section 
containing the waveguides N + 1, TV + 2, ..., is differ- 
ent. As we demonstrate below, the variation of the back- 
ground coupling enables selective manipulation of Fano 
modes and regular exponentially localized modes. 

We model the evolution of the optical mode amplitudes 
along the propagation direction using a coupled-mode ap- 
proach with nearest-neighbor coupling [25] , 

i-^- + CWi+i + t/>j-i) = 0, 1 < j < N 
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^ + C 6 (^+i+^-i)=0, 3>N + 1. 

where ipj describes the optical mode amplitude in the 
j-th waveguide. We put ifj = due to the array ter- 
mination. In this formulation, the spectrum is centered 
around the propagation constant of an isolated waveg- 
uide. 

We first consider the case of a homogeneous coupling 
in the waveguide array, Cb = (7. Then, the eigenvalues j3 
in the semi-infinite array for a continuous band in the re- 
gion j3 G [—2(7,2(7]. By virtue of the dispersion relation 
/3/C = 2 cos(ft), each supermode can be constructed from 
plane waves exp(mn) with wavenumber n. The waveg- 
uide F is a defect to the ID chain and plane waves in- 
cident on the triangle of waveguides exhibit a transmis- 
sion which is dependent on the value of k. It is con- 
venient to represent the field inside the resonator region 
(n = 1, . . . , N) in terms of incident and reflected waves as 
ip n = exp(i/3z)[exp(i(n — N)n) + i?exp(— i(n — N)k)] and 
in the background region (n = iV+l,iV + 2,...) the field 
will be the transmitted wave i/j n = ex.p(if3z)T ex.-p[i(n — 
N - 1)k]. One finds from Eq. ^ with C b = C that the 
transmission and reflection coefficients are 



(2) 



T = ie iK [C 2 sin(3«/2) + (V 2 - C 2 ) sin(«/2)] D~\ 
R = e iK cos(K/2)V 2 D-\ 

where D = C 2 [zsin(3^/2) + cos(«/2)] + (V 2 - 
C 2 ) exp(m/2). We note that the transmission vanishes at 
certain value of k = ftp, which corresponds to a Fano res- 
onance [25]. At resonance, we have R(k f ) = exp(iftp), 
i.e. Kp defines the phase of the reflection coefficient. In 
order to use the Fano resonance to trap light in the res- 
onator region, it is necessary to also satisfy the Fabry- 



Perot condition for one round-trip of the plane wave be- 
tween the surface of the chain and the AW, which is for- 
mulated as 

2(7V - 1)k f + n F + 2k f = k f (2N + 1) = 27rra, (3) 

where m is an integer. In this expression, the first term 
corresponds to phase accumulation for the waves propa- 
gation between the waveguides n— 1, . . . , N, the second 
term is the reflection phase from the Fano defect, and 
the third term is the reflection phase at the boundary 
n = l. We find that the Fano and Fabry- Perot reso- 
nance conditions can be satisfied simultaneously when 
V — C. In this case k F = 27r/3, and the Fabry- 
Perot resonances occur for m = 1,3,9,..., correspond- 
ing to the number of waveguides in the resonator region 
A/" = l + (3/2)(m-l) = l,4,7,... 

In the following, the case N = 4 is studied numer- 
ically and experimentally. For the experimental study, 
the waveguide array has to be finite, and we consider an 
array with all in all J = 49 waveguides (48 waveguides in 
the ID chain plus the AW). We numerically calculate the 
spectrum of eigenmodes, and present it in Fig. [2|a). We 
find extended eigenstates with eigenvalues f3 G [—2(7, 2(7] 
representing a quasi-continuum. Furthermore, as indi- 
cated by arrows in the figure we identify two localized 
states, which arise from the coupling of the AW with the 
ID chain. The state indicated by a red arrow is a defect 
state residing inside the gap above the band of eigenval- 
ues and it is thus exponentially localized (/? = 2.38(7). 
The second state, indicated by the blue arrow, occurs at 
the eigenvalue corresponding to the Fano resonance, at 
it has a fundamentally different field structure with the 
zero amplitude in the background, i.e., this state strictly 
exists only between the edge and the AW. This prop- 
erty is referred to as "compactness". Furthermore, the 
propagation constant of this Fano-compact state (FCS) 
resides inside the allowed propagation band. The Fano 
resonance does not break up the band of eigenvalues of 
the unperturbed ID chain which makes the FCS a bound 
state in the quasi-continuum of the ID chain. The van- 
ishing of the electric field in waveguides N = 5, 6, 7, . . . 
does not allow any interaction of the FCS with the resid- 
ual lattice. Therefore the state is completely insensitive 
to the lattice structure in the background. Hence, the 
FCS does not change if Cb or the number of waveguides 
in the background is varied. Figure |2|b) demonstrates 
the change of the FCS (left) and the defect state (right) 
with increasing C&. By an increase of Cb > C the prop- 
agation band broadens until the exponentially localized 
state is absorbed into the band, whereas the FCS re- 
mains unchanged at e F = — 1. The critical ratio of the 
coupling constants for which the exponentially localized 
state ceases to exist is C£ r jC = 1.36. 

A light beam coupled to the waveguide array will excite 
a superposition of different modes. The amplitudes of 
eigenstates will be proportional to the overlap of their 



2.38C 

(a) 2C 

CO. 

3 

to 

I o 

£= 
.O 



Q. 

o 



-2C 



f^" 



F 



OOOOOI 



3 

o 



2 

o 



electric field amplitude (arb. u.) 



F 



o o o 6 o o 



3 

o 



2 



compact state 



defect state 




waveguide no, 



waveguide no. 



FIG. 2. (Color online) (a) Eigenvalue spectrum of a struc- 
ture with 49 homogeneously coupled waveguides (Cb = V = 
C). The amplitude profiles of the exponentially localized de- 
fect state (upper row) and the FCS (lower row) are shown on 
the right. (b,c) The profile of (b) FCA and (c) defect state as 
functions of C&. The limit where the exponentially localized 
state ceases to exist is C£ r = 1.36C. For Cb > C^ r the defect 
state is absorbed in the quasi-continuum and hence spreads 
over the entire array. 



spatial profiles with the input field distribution. For a 
single- waveguide excitation of the first site in the above 
structure, the amount of intensity dedicated to the FCS 
is 25%. The defect state is only excited with a fraction 
of 0.37% of the input intensity when C& = C, and this 
decreases further when Cb > C. 

For our experiments, we fabricated several waveguid- 
ing structures in fused silica using the femtosecond laser 
direct- writing approach [26]. The samples are 10cm long, 
which is the maximum length feasible with our current 
fabrication technology. The experimental investigations 
now deal with two separate issues. First, observing the 
Fano resonance in an array with homogeneous coupling 
in the entire structure. Second, proving the compact- 
ness of the FCS by varying the coupling constant Cb in 
the background. In order to directly observe the light 
evolution inside the structures we employ a fluorescence 
microscopy technique [27 . Among the possible single- 
waveguide excitations, the excitation of the first waveg- 
uide at the edge of the structure has the largest overlap 
with the FCS compared to the overlap with the defect 
state. In our experiments light at A = 633nm is launched 
into this waveguide using fiber butt coupling. In a first 
set of experiments the spacing between neighboring wave- 
guides of the chain is fixed to 13/im, such that Cb = C. 
This results in a coupling of about C = 0.2mm" 1 . Fig- 




FIG. 3. (Color online) (a) Microscope image of the front 
facet of the processed fused silica sample serving as the ex- 
perimental realization of the structure, (b) Corresponding 
fluorescence microscopy measurement and the near-field im- 
age of the output facet, (c) Computed intensity evolution for 
20cm of propagation. The fluorescence images are normalized 
to their respective peak value. 



ure[3ja) shows the microscope image of the front facet of 
the fabricated array. Due to the strongly elliptical shape 
of the waveguides, resulting from the fabrication process, 
the coupling strength depends not only on the separation 
of the guides but also on their orientation with respect 
to the coupling direction [28]. The vertical offset (h) of 
the AW with respect to the chain was carefully tuned to 
match the condition V = C. Equal coupling was achieved 
for h = (V3d/2 + 2.6)/xm. 

The experimental data is shown in Fig. |3|b). When 
launching light into the edge waveguide of the structure, 
several eigenstates of the system are excited, including 
the FCS. Due to the reflection at the defect caused by 
the AW, the light carried by the extended eigenstates is 
transmitted gradually into the background. After 10cm 
propagation the intensity pattern in the resonator region 
is not perfectly reduced to the intensity pattern of the 
FCS. The transient behavior within the first five wave- 
guides still includes intensity oscillations in each waveg- 
uide. However, the output pattern already exhibits the 
characteristic "dark" third waveguide [Fig. [3lb), lower 
part] identifying the FCS. The simulations, shown in 
Fig. [3|c), confirm the measured behavior. Even after 
20cm propagation, the FCS is not established perfectly 
but the data in Fig. [3] is sufficient to confirm the existence 
of the FCS. In the next experiment, we aim to demon- 
strate the "compactness" of the FCS. To this end, in a 
second sample we reduced the spacing in the background 




FIG. 4. (Color online) Light evolution in a structure, 

where the coupling in the background region is increased to 
Cb/C =1.75. (a) Fluorescence microscopy measurement and 
the near- field image of the output facet, (b) Computed inten- 
sity evolution for 20cm of propagation. Only the fluorescence 
images are normalized to their respective peak value. 



region to 11/im which corresponds to C& ~ 0.35mm -1 , 
whereas in the region close to the surface 13/im spac- 
ing is reproduced. Therefore, the Fano resonance is ex- 
pected at a similar height h of the AW, compared to 
the first structure. Indeed, V = C was here achieved for 
h = (v // 3^/2+2.2)/im. The experimental data in Fig. Ufa) 
shows that an intensity distribution in the first five wave- 
guides evolves to a final state which is very similar to that 
observed in the first experiment, c.f. Fig.J3^a). Our simu- 
lations of the light evolution in this structure [Fig. [4^b)] 
confirm this trend. This is the experimental confirma- 
tion that the FCS is independent of the topology of the 
background region of the ID waveguide chain and, thus, 
the FCS is compact. Nevertheless, the transient behavior 
changes slightly. This is reasonable since the occupation 
of the eigenmodes has changed and except for the FCS, 
the profiles of eigenmodes are sensitive to the value of 

c b . 

In this work, we have theoretically predicted and ex- 
perimentally demonstrated the existence of a compact 
bound state in the continuum of one-dimensional waveg- 
uide array. This state is confined to a finite region of the 
array and exhibits no exponentially decaying tails. We 
found the condition when all localized states other then 
the compact state cease to exist. 
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Methods 

Our samples are fabricated in bulk fused silica wafers us- 
ing the femtosecond laser direct- write approach [26] em- 
ploying a Ti: Sapphire Mira/RegA laser system (Coherent 
Inc.) operating at a wavelength of 800 nm, a repetition 
rate of 100 kHz and a pulse length of 170 fs. The light 
is focused inside the sample by a 20 x microscope objec- 
tive (NA=0.35), and by continuously moving the sample 
using a high-precision positioning system (Aerotech Inc.) 
the waveguides are created by the induced refractive in- 
dex increase. For the fabrication of our samples, the pulse 
energy was adjusted to 320 nJ and the writing velocity 
was set to 1.5 mm/s. 
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